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Abstract 

In the present work a mathematical model aimed at evaluating the dynamics of crowds is 
derived and discussed. In particular, on the basis of some phenomenological considerations 
focused on pedestrians' individual behavior, a model dealing with the collective evolution of 
the crowd, formalized in terms of a conservation law for a crowding measure, is obtained. To 
suit engineering needs, modeling efforts are made toward the simulation of real world crowd 
events specifically happening along pedestrian walkways. The response of the model as well 
as its tunability are studied in some computational domains inspired by existing footbridges. 
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1 Introduction 

The study of the dynamics of crowds defines a wide research sector that, in recent times, is know- 
ing an increasing expansion fostered by different interests and disciplines. In the authors' opinion, 
throughout the years, the multiple facets of the field have been approached almost independently 
by different scientific communities, such as those of Applied Mathematics, Physics, Biomechan- 
ics, Transportation and Civil Engineering [29]. The present study develops from a collaborative 
multidisciplinary approach between applied mathematicians and structural engineers and aims at 
crossing and linking such facets. First, a mathematical model to simulate the evolution of crowds 
is introduced and inquired into; second, modeling efforts are made toward the simulation of real 
world crowd events specifically happening along footbridges. 

In the proposed model, the crowd evolution is primarily approached considering the phe- 
nomenological perspective of an individual pedestrian. In particular, an evolution equation for 
agents' positions is derived. Such an equation takes into account the active attitude of pedestri- 
ans, namely the fact that their motion is not passively driven by force fields, as it happens with 
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inert particles. Rather, pedestrian positions evolve according to personal desires and external stim- 
uli (see e.g., [Toj[T2j[T6] ; the reader can refer to |2,14 for a general picture on crowd modeling). 
However, since macroscopic or statistical information is more significant in the engineering design 
and evaluation of the facility performances, a formal procedure is introduced to derive an evolu- 
tion equation for pedestrian collectivity. Such a procedure, in the spirit of Reynolds Transport 
Theorem, generates an evolution equation for a mesoscopic measure of the pedestrian crowding, 
in the form of a conservation law. Moreover, in order to meet engineering codified approaches, a 
statistical interpretation of the model output is further proposed. 

To deduce the model, completeness as well as minimality are assumed as conceptual guidelines. 
In particular, two kinds of interactions are retained for determining the dynamics: the interac- 
tions among pedestrians and those between pedestrians and domain walls. Doing so allows one 
to obtain simple equations depending on a limited amount of tunable free parameters. Further- 
more, looking from the engineering practice perspective, ways of imposing boundary conditions, 
handling articulated domain geometries, and discretizing the equations of motion are proposed to 
complement the mathematical model. Finally, the response of the model as well as its tunability 
are studied in some computational domains inspired by existing footbridges. 

The paper is organized in four more sections that follow the conceptual approach described 
above. In particular, in Sect. [2] the phenomenological model is proposed and the mathematical 
model is therefrom deduced; in Sect. [3] further elements required to simulate real world crowd 
events along footbridges are discussed; in Sect. [4] results from simulations of crowd events in four 
computational domains inspired by real world as well as the tunability of the model are discussed. 
Finally, in Sect. [5] conclusions are outlined. 



2 Crowd modeling 

The aim of this section is to introduce and discuss the crowd model considered throughout the 
paper. A phenomenological model describing the dynamics of pedestrians in a moving crowd is 
firstly obtained (Sect. |2.1[ ). On such a basis, a mathematical model suitable to provide data on 



the collective evolution of the crowd is derived (Sect. 2.2). Finally, a statistical interpretation 
of the latter is discussed (Sect. |2.3[ ). This interpretation is aimed at providing statistical data 
on crowds, in the direction of the current engineering practice, which, more and more frequently, 
require probabilistic characterization of events (see e.g., |1 in Civil Engineering). 

In the modeling approach we propose, individual pedestrians actively determine their velocity 
on the basis of interactions they have with the collectivity around them, being the latter observed 
at a mesoscopic-\ike scale. Indeed, in this approach, the individual and the collective points of 
view coexist and interact, determining the overall system dynamics. For the sake of simplicity, 
the pedestrian-collectivity interaction mechanism is treated considering two phenomenological 
aspects: the perception a walker has of her close neighbors and her consequent reaction. From 
the modeling point of view, the coupling of these two aspects is done via a velocity term called 
interaction velocity. Such a term is calculated on the basis of a spatial measure^ namely the 
mesoscopic variable involved, which models what one can perceive. In the model, the interaction 
velocity plays the role of an additive perturbation of a second velocity term, the desired velocity, 
which, instead, models the velocity one would keep in the absence of other nearby agents. The 
sum of desired velocity and interaction velocity gives rise to pedestrian total velocity. 



2.1 Phenomenological modeling of pedestrians dynamics 

Let us imagine a crowd walking in a domain D CW^. Usually d = 2, but the approach we present 
is sufficiently general to allow one to model possibly different systems of active interacting agents 
in different space dimensions. For instance, vehicular traffic with d — 1 or swarms with d = 3. 
Therefore, we keep the dimension d generic. Throughout the current section we set D = R^; 
nonetheless, in Sect. [3] this requirement will be relaxed allowing bounded domains, closer to 
engineering needs, to be considered. 
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At time t > 0, the point of view of individuals is treated in terms of the spatial position Xt ^ D 
of a generic pedestrian, henceforth referred to as test pedestrian. On the other hand, the point of 
view of collectivity is expressed in terms of a time evolving positive Radon measure Vf Such a 
measure, supported on quantifies the pedestrian mass in any Borel measurable set E G B{D)^ 
being B{D) the Borel a-algebra of i.e. 

yt{E) = measure of the crowding of E at time t. ^ ' 

From now on, we will simply write E C D to denote at once a subset of D which is also measurable, 
i.e., belonging to B{D). 

Owing to Lebesgue Decomposition Theorem, the measure admits a unique decomposition 
as a sum of an absolutely continuous part and a singular part with respect to a given measure, 
which, in our case, is the d-dimensional Lebesgue measure in W^. Moreover, by Radon-Nykodim 
theorem, the absolutely continuous part can be expressed by means of a non negative density 
function pt G L^{D). On the other hand, the singular component can be further written as an 
atomic part (i.e., a sum of Dirac masses denoted by u^) plus a Cantor part (see e.g., |22 ). In 
view of modeling human perception, we are essentially interested in the absolutely continuous and 
atomic parts, hence we set: 

d^t = ptdx -\- dv^ . 

In particular, the density pt is intended to model a gross, coarse grained, perception of individuals, 
while the atoms constituting are aimed at reproducing a more granular and sharper perception, 
which is precise and spatially localized. 

As previously mentioned, pedestrians are assumed to have a desired walk velocity, which is 
kept in absence of other nearby people. This velocity is expressed through a vector field 

which is evaluated at the agent's position. On the other hand, to obtain the interaction velocity^ 
a mapping 

k'.DxD ^M."^ 

is considered, so that /c(Xt, y) models the reaction that the test pedestrian in Xt has to a pedestrian 
perceived in y. As a pedestrian generally interacts just with nearby agents, on the basis of what 
she perceives^ the interaction velocity is recovered first by restricting /c to a region S{Xt) G B{D) 
around X^, and then by weighting the pairwise interaction expressed by k with the perceived mass 
vt, i.e. 

v,[vt]{Xt)= I k{Xt,y)dut{y)= f Is^x^K^u v) dvt{y). (2) 
JS{Xt) Jd^ V " 

■.=KiXt,y) 

where ls{Xt) denotes the characteristic function of the set S{Xt). The latter is intended to model 
to the so-called sensory region of the test pedestrian (see Fig. [l]and, e.g., [lO]). Coherently with 
human physiology, such a region is expected to be bounded and possibly anisotropic with respect 
to the pedestrian direction. Since the sensory region is bounded, agents cannot voluntarily produce 
the collective trends that one can observe by looking at their global distribution because they 
have a limited perception of the group they are part of. 

In order to obtain the agent's total velocity, the interaction velocity is added to the desired 
velocity, giving 

dX f 

= VdiXt) + Vi[vt]{Xt) = VdiXt) + KiXt, y) dvt{y), (3) 

which details how Vi generates deviations from the trajectory indicated by Vd- 

The model deduced thus far features two state variables, Xt and Vti although just one evolution 
equation ^ has been provided. A second equation can be obtained via the following argument. In 
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Figure 1: Sensory region of a given pedestrian - pedestrians perceived in terms of their distributed 
(checker texture) or locahzed (arrow) mass 



order to be consistent with its definition of pedestrian mass, Ut must be a material quantity with 
respect to pedestrians and their movement. In other words, the initial perceived mass distribution 
uq has to be transported by pedestrians themselves. This fact is expressed by the equation 

Ut=Xt#uo, (4) 

which states that Vt is the push forward through Xf of the initial perceived pedestrian mass 
distribution uq. Equation Q means that for every measurable subset E of the relation 

holds. 



2.2 Mathematical model 

In this section, we derive a self consistent mathematical model from the coupled evolution equations 
([3| and Q. To this end, we compute a formal expression of the time derivative of the measure 
z^t, thereby obtaining a conservation law which expresses its evolution in time and space. 

Let us consider Ut as a distribution, i.e., for every smooth test function tjj G C^{D) compactly 
supported in D we consider the action 

JD 

In view of Eq. Q, the duality (?/^, Vt) satisfies 

{^.yt)= I ^d{Xt#iyo)= [ ^{Xt)diyo^ (5) 
Jd j d 

hence, taking a formal time derivative (very much in the spirit of Reynolds Transport Theorem), 
we recover 

From Eq. ([3| we further have 

VV(Xt) ■'^dvo = (^^{Xt) ■ [vd{Xt) + K{Xt, y) diyt{y)^^ di^o, 
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thus 

vt) = {^d + J^^(-^ y)dMy?j^ dvt. (6) 

We say that the measure Vt satisfies in the weak (i.e. distributional) sense the conservation la^^ 

^+V-{vt{va + Vi[vt]))=Q (7) 

provided Eq. ([g]) holds for all test functions G C^{D) and for almost every (a.e.) time t in a 
given time interval of interest. 

Under reasonable assumptions on the velocity field, one can have existence of solutions to the 
Cauchy problem obtained by complementing Eq. ^ with an initial condition 

uo = V. (8) 

The interested reader is referred to for technical details. 

The great generality of the concept of measure has been transferred into Eq. ([7|). For instance, 
if Vt is the counting measure, the obtained system is formally equivalent to a discrete dynamical 
system. On the other hand, Vt can be absolutely continuous with respect to the Lebesgue measure 
on D and thus ^ becomes a usual conservation law for a density. As already mentioned, Vt can 
also be hybrid in nature, having a discrete and a continuous component at the same time. This 
allows one to embed the dynamical effects of the microscopic granularity typical of crowds also in 
a continuous description given in terms of a density f8|. 



2.3 Statistical interpretation of 

The model deduced thus far is self-consistent and can be considered on its own. However, a further 
statistical reading can be attached to it. Such a reading has a twofold aim. Firstly, it allows one 
to better inquire into the real scale at which we are observing the crowd, i.e., the scale which 
the measure Ut refers to. Secondly, statistic-like information is required in the current engineering 
practice as previously mentioned. 

Let us suppose that the considered crowd is constituted by N indistinguishable agents. Given 
an abstract sample space {ft, T ^ P), we consider the N random variables Xl^Xl^ . . . , referring 
to their position, i.e., 

Xl'(Sl,T,^)^D, VjG{l, 2...,7V} 
Xl = position of agent j at time t. 

Of course, the X^'s, being positions of pedestrians, evolve according to ([3|. 

Let us further assume that, at time t = 0, the positions Xq identically distribute in D according 
to the normalized initial perceived mass measure, i.e., 

n*^ = (9) 

:=/x 

that is 

P(X^' eE) =fl{E), \fEcD. 

Equation (|9| can be motivated, at least intuitively, on the basis of the idea that the probability 
of finding a pedestrian is higher where more pedestrians are present. 

Because X^ solves Eq. ([3|, the following identity between measurable mappings from Q into 
D holds: 

X/(cj) = {Xt o X^q){uj), for a.e. uj e (10) 

(see also Fig. [2]). This equation details the real nature of the stochasticity of the X/'s, which is just 
a legacy of the stochasticity in the initial positions Xq. In fact, the evolution of the uncertainty 

■"^The term Vi has been defined in Eq. ([2J 
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Figure 2: Conceptual representation of the measurable maps Xq and Xl 

on the X^'s at any time is deterministically operated by the flow map Xf. In particular, the 
probability measure Xl #P satisfles 

Xi#W = {Xt o X^')#P = Xt#{Xim = Xt#Ji = ^Xt#v = ^vt. 

Hence, the measure jit j^^t is exactly the probability law referring to the position of the j-th 
pedestrian at time t, i.e. 

P(X/ eE) = fit{E), WE CD. 

In other words, the X^'s, "dragged" by X^, transport the uncertainty associated to their initial 
condition. 

Accepting this statistical interpretation, one can count the total amount of pedestrians in a 
given set E C D at at time t through the variable 

N 

which is random because it is a function of the random variables X^. Its expected value can be 
easily computed: 

AT 

nyEA := ^E[l^(X,^')] = N^,{E) = Vt{E). (11) 

The quantity Vti initially introduced as a mesoscopic quantity referring to collectivity, can clearly 
have a macroscopic interpretation as an expression of the average pedestrian mass. On the other 
hand, a linear relation exists between Vt and /i^, and the latter embeds a genuinely mesoscopic-like 
information. Actually, this scaling is linear because we are considering flrst order models, i.e., 
models which do not invoke explicitly any concept of acceleration. The link just found can have 
some flavor of "equipartition of energy", namely all pedestrians contribute the same way to the 
macroscopic mass as they are indistinguishable. 

3 Toward the simulation of real world crowd events 

The aim of this section is to obtain a tool to analyze crowd events in real domains on the basis of 
the modeling framework previously deduced. From now on, two working assumptions are made: 
the measures involved are assumed to be absolutely continuous with respect to the Lebesgue 
measure, hence any granular effect is neglected, and the dimension of the working domain is flxed 
io d = 2. 

First, to obtain a self-consistent tool a functional expression must be given to the interaction 
kernel deflned in Eq. In order to meet engineering needs, such an expression is expected to 
be synthetic, although obtained on a phenomenological basis. 
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Figure 3: Sensory region S^{Xt) of the test pedestrian in Xt 

Second, the current study considers problems concerning crowd flows in built environment. 
Such flows usually develop in very large phenomenological environments, in the sense that pedes- 
trian dynamics are affected by phenomena which originate far from the spatial region of engineering 
interest]^ Thus, on the one hand, the mathematical domain is a restriction of the phenomenolog- 
ical environment and inflow conditions need to be prescribed on some conventional boundaries. 
On the other hand, the environment is bounded by built perimeters having articulated geometries 
in plan (e.g., building facades, obstacles, parapets). Hence, pedestrian behavior must be sensi- 
tized to walls and the discretization procedure must be able to handle such possibly complicated 
geometries. In the following, our work will be mostly inspired by problems arising when consid- 
ering footbridges. In particular, the computational domains considered own to the class of the 
elongated domains and the considered inflow conditions may mimic queue processes happening on 
footbridges [TT]. 

In summary, modeling solutions have to be given so as to address the following aspects: i) in- 
teractions among individuals; ii) limited spatial domains D CM?; iii) pedestrian behavior at walls; 
iv) numerical implementation in complicated domain geometries; v) crowd inflow. The current 
section details the modeling solutions for these aspects. 



3.1 Modeling pedestrian interactions 

The interaction velocity Vi introduced in Sect. Q is the modeling artifact to account for the 
interaction among pedestrians. The kernel K = K{Xt^ y), in particular, is in charge of modeling 
its intangible psychological aspects. It is worth remarking that, in order to make model tuning 
affordable, K should depend on a limited amount of free parameters. 

In the current section, an expression of K based on three positive parameters is proposed. In 
Sect. [4j the parameters involved will be tuned and used to perform simulations. 

Commonly, in normal crowd avoidance, pedestrian interactions are repulsive and limited to a 
frontal region (being, thus, anisotropic). Therefore a minimal expression of K can be: 

KiXt,y) = -c^ls^^^x,)iy), (12) 

where c is a positive constant (which determines the repulsiveness of the interaction) and r = 
\y — Xt\ is the inter-pedestrian distance whose direction is = \yZxl\ • 



The interaction described in Eq. (12) is conflned to the sensory region /S(X^), which we model 
as a circular sector of radius center X^, and angular semi-amplitude < a < 7r/2 around 
V(i{Xt). Formally, the set S{Xt) can be written as 

S{Xt) = S%{Xt) = {yeD :r<R, Vd{Xt) • > \vd{Xt)\ cos a} , 

see Fig. |3] 

^For instance, one can consider spectators leaving a stad ium at the end of a sport event and then crossing 
a footbridge which connects the stadium with the car park |ll]. A structural/transportation engineer may be 
interested in simulating the pedestrian flow over the footbridge, omitting a detailed description of the grandstand 
evacuation process. 
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Figure 4: (a) Rectangular domain D - (b) Detail of D with the desired velocity vector field, c) 
Iso-lines of angle 7 



3.2 Modeling the effects of walls 

The model considered thus far operates on the whole unbounded physical space M?. However, 
domains relevant in applications commonly feature built perimeters, in general walls^ which agents 
cannot cross. Therefore, proper behavioral rules describing the agents' reaction to walls should be 
introduced. 

The rules proposed in the current section are obtained on the basis of a pure physical intuition, 
supported nonetheless by the evidence that both desired and total velocities should not allow for 
wall penetration. In particular, constructing suitable desired velocity fields in presence of walls 
for generic bounded domains is a nontrivial task, thus in the next section elongated domains will 
be specifically considered in view of the application to footbridges. 



3.2.1 Modeling the desired velocity in presence of walls 

In the current modeling framework, the agent desired velocity field Vd, which ideally drives pedes- 
trian motion in the domain, is supposed to be known a priori. Therefore, it should be constructed 
out of the knowledge of the geometry of the domain. 

Let us consider a generic wall bounding the domain. We denote by n the outward normal unit 
vector (see Fig. [4jb)). Because of the impenetrableness constraint, the following compatibility 
condition, to be understood as a basic design guideline for can be introduced: 

Vd-n<0. (13) 

When general domains are considered, constructing a field Vd which is phenomenologically 



acceptable and satisfies Eq. (13) is a nontrivial task, which has been considered both in crowd 
modeling literature [T7j[T9] and, more generally, when path planning is concerned (e.g., robot 
motion planning f^). In the following, a method to build such a field in case of elongated domains 
is proposed. These domains are characterized by a longitudinal dimension crossed by pedestrians 
(length) much larger than the transversal one (chord), the latter possibly slowly varying along the 
former. They may model, for instance, pedestrian footbridges, sidewalks, platforms and so on. 
Moreover, as in [T7||T9], we wish to consider velocity fields that are (normalized) potential fields, 
i.e., such that Vd = — Vi^/|Vi^| for some potential function u that has to be determined. 

Let us initially consider the simplest domain geometry, where the problem can be easily set 
and modeling considerations appear more intuitive. Let D be a rectangular domain of length L, 
chord and aspect ratio B = B / L = ^ (Fig. [ij^a)). To find we assume: 

• that pedestrians fiow from left to right in D. Therefore, the field —Vu has to be rightward. 
Moreover, a unit potential difference across L is assumed; 




• that pedestrians aim at not scraping against lateral boundaries. Hence, —\/u must be 
directed inwards at the ends of the chord and longitudinally at mid-chord. In other words, 
being Cx the longitudinal direction, the angle 7 = cos~^ {{vd ■ ex) / \ vd\) is expected to decrease 
monotonically to zero when approaching the mid-chord. We parametrize the convergence 
rate of Vd in terms of its slope at walls: 

tan6> = tan7|r,, (14) 

see Fig. ^h). We stress that, since inter-pedestrian interactions are repulsive, if pedestrian- 
wall repulsion is not taken into account then agent agglomerations are likely to appear in 
the proximity of the lateral boundaries Tg. 

A minimal potential complying with the assumptions above is: 

u{x, y) = -x^qy'^, (15) 

where the coordinates x and y are scaled with respect to the span L and moreover q = tan6>/5. 
Hence, u defines the potential field —Vu = (1, —2qy)^ which, up to normalization, generates the 
desired velocity field Vd = (1, —2qy)/^s/l + Aq^y'^. 

In order to deal with more general low aspect ratio domains, we primarily observe that Eq. 



( 15 ) solves the following Poisson problem: 




-2q in D 

= tanl9| onT, (16) 
= —Xi + qy on Ti {i G {m, out})^ 

which can be imposed on more general domains than the rectangle. The term is intended to 
parametrize the convergence rate of Vd toward the longitudinal direction at mid-chord. Therefore, 
the factor h/B^ where h = b{x) is the (spatially variable) chord amplitude, is introduced. Factually, 
in rectangular domains, where b = the Neumann boundary condition on Tg turns out to be 



exactly condition ( 14 ) 



The potential defined in Eq. (16) is subharmonic. In the literature, analogous velocity fields 
have been usually built by using harmonic functions [7|[T9j[2l], relying on the so called min-max 
principle. Such a principle ensures that any local maximum or minimum of i.e., a stationary 
point of the gradient field, lies on the boundary of D Therefore, provided proper Dirichlet 
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conditions are imposed (e.g., = on exits, u = 1 anywhere else on the boundary), velocity 
fields driving pedestrians to a given exit can be obtained. On the other hand, no control on the 
behavior of —\/u at walls is allowed, hence phenomenologically unsatisfactory behaviors can arise. 
Nevertheless, subharmonic functions do not satisfy a minimum principle, hence the absence of 



internal local minima cannot be guaranteed anymore. Aside from this, their use through Eq. (16) 
allows phenomenologically consistent fields |15 to be obtained. 

The model in Eq. (16) with ^ = 5° is applied to some real world footbridges having different 
walkway shapes (Fig. ^5f: bottleneck shape (Chiaves footbridge, |23 ), curved shape (Liberty 
footbridge, |3 ), shifted shape (Coimbra footbridge, |6 ). As far as the longitudinal direction is 
concerned (Fig. [5|2)), the obtained fields are constant and a unit potential difference across L is 
achieved. For what concerns the chord-wise direction, the desired walking angle 7 is sensitized to 
the local sidewall geometrical inclination, for the sake of clarity. To do so, we introduce the angle 
P := 7 — a^, where extends the geometric angle between the boundary and e^^. In particular, it 
is expressed in terms of the function ai{y) = \ai{y — y2) / {yi — ^2) — OL2{y — Hi) / {y2 — yi)\^ where 
aj {j = 1, 2) are defined in Fig. |5|2a). In Fig. [5|3) the angle /3 is plotted along the sections cl 
{x = — 0.35L), c2 {x = 0, mid-span), c3 {x = 0.35L) (Fig. |5|2)). It is worth pointing out that 
the obtained (3 chord- wise profiles are close to the expected trend (gray dotted lines. Fig. [5|3)) 
when the walkway section is almost constant; nonetheless, the /3 profile coherently departs from 
this trend when significant geometry variations take place (e.g. |5|a)-c2, [5|c)-c2). 

3.2.2 Effects of walls on the total velocity 

The constraint of not penetrating walls refiects also on the conditions to be enforced on the total 
pedestrian velocity at the boundaries. The following conditions should then be satisfied: 

^•n={vd^Vi)-n<0. (17) 

Two complying conditions, indeed referring to antipodal behaviors, can be found: 
1. Pedestrians may scrape against walls, i.e., 

= {vd + Vi) - (n • {vd + Vi) n) [{vd + ^i) • n > 0] ; 

walls 



dt 

2. Pedestrians must stop at walls, i.e., 

= {vd^Vi)[{vd^Vi)-n<0]. (18) 



dt 



walls 



Equation (18) can be thought as an equivalent of the no-slip condition in fiuid mechanics. In the 



simulations presented in Sect. [4] condition [T] is used. 



3.3 Simulation of crowd events in articulated domains 

The current section deals with the numerical discretization of Eq. ^ for producing numerical 
simulations of crowd events in real, possibly articulated, domains. Spatial discretizations, done 
via unstructured triangular meshes, are considered. 

Problem can be put in discrete form by using the ad hoc scheme introduced in |19 . 

Operatively, discretized equations are obtained after a twofold discretization in time and space. 
First, a first order explicit-in-time discretization is needed. Let [0, T] be a time interval of interest, 
which we evenly partition by the instants tn = nAt, so that n ranges from to a value M such 
that T = MAt. In this setting, we recursively generate a countable family of measures {i^n}^=o 
approximating the i^tj^- To this end, let us introduce the the discrete-in-time fiow map: 

Xn{x) =X^Wn{x)At, (19) 
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Figure 6: Conceptual sketch of the numerical scheme in action 



where Wn{x) := Vd{x) + Vil^njix). By using Eq. (19), one can define the recursion 



\ i>o = 

which produces the desired time approximation: 

z>n ^ J^nAt, n = 1, 2, . . . , M. 

As u is absolutely continuous, if the discrete flow map meets suitable regularity conditions then i>n 
is absolutely continuous as well for all n (for a detailed discussion on the requirements the reader 
can refer to |19 ). 

After the time discretization has been established, we consider a finite-volume-type partition 
of D. Let {Ej.}^^-^ be a grid of Q measurable elements of centroids and let pn be the density 
of i-e., di>n = Pndx. We approximate pn by means of a piecewise constant function pn given 
by: 

Q 

Pn{x) = ^p11e^{x), 

k=l 

where p^ is a characteristic value of pn when restricted to the element Ej. (e.g., p^ = pn{x^))' We 
finally call z>n the measure whose density is pn, i.e., dOn = pn dx. 



A piecewise constant space approximation of the flow map (19) can be further considered: 

^n{x) = X -\- Wn{x)At^ 

where the piecewise version of the velocity Wn{x) is used, i.e.. 



W. 

k=l 

As the mesh is fixed, i.e., it does not get deformed in time, the recursive relation 

Z>n+1 = ^n#Z>n, n > 1 

is further approximated by 

Q 



^0 _ v-^fe l-E^gnX,^(Efc)| o 

— 2-^ Pn \ 5 ^ — J-, ^, • • • , V, 



(20) 



where, for the sake of conciseness, we have denoted by | • | the Lebesgue measure . If the spa- 
tiotemporal grid is properly refined (viz. under a suitable relationship between the characteristic 
size of the elements and the time step At), the scheme converges to the weak solution of Problem 
0-(|8| (for technical details see [26 ). 

This numerical scheme can be ideally used with any type of grid, independently of the element 
shape. From the strict implementation point of view, one needs to compute rigid movements of the 
elements Ek and evaluate their intersections. For instance, in |8|n!9] orthogonal grids formed by 
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Figure 7: Pedestrian entering region X 



square- shaped elements have been used. Nonetheless, when articulated domains are considered, 
triangular meshes are often used. Pairs of triangles exhibit a wide selection of (topologically) 
different ways of intersecting (see Fig. [g] and e.g., |24 ), which makes the evaluation of ([20| 
expensive. It is worth pointing out that the problem of intersecting triangles (and convex polygon 
in general) is typical in computational geometry and computer graphics and can be solved very 
efficiently. Algorithms having linear complexity with respect to the number of edges of the polygons 
have been conceived (see e.g., 1 18,27 ), which allow computational times to be contained. Domains 
considered in simulations in Sect. [4] (see Fig. |5| have been discretized by using triangular meshes. 



3.4 Handling pedestrian entering 

Often, in applications, the spatial domain considered is just a subset of a larger phenomenological 
domain in which the whole crowd event develops. Therefore, assuming that at the initial time 
all involved pedestrians are already inside the domain can be over-restrictive. This is especially 
true when footbridges are considered. On the contrary, this is not the case when dealing, for 
instance, with evacuation problems [25 . In the current section, an approach allowing for pedestrian 
entrance, conceptually consistent with the model previously developed, is proposed. Specifically, 
arrival processes mimicking queues will be considered. 

As in Fig. [t) let us consider a region X C adjacent to D which models a buffer through which 
pedestrians coming from an external (i.e., not explicitly modeled) part of the phenomenological 
domain enter D. In other words, X can be thought of as a lobby or a waiting room. Consistently, 
we assume that it has a non-zero Lebesgue measure, i.e. \X\ > 0. From now on, X will be referred 
to as the entering region. 

For our purposes, let us imagine the perceived mass measure of X as a bulk, super-macroscopic, 
variable of the system, which we henceforth denote by It := i^t{^)- At time t, we denote by St 
the number of pedestrians involved in the crowd event who are waiting to appear in X. They are 
ideally stored in a 0-dimensional reservoir S. Once they are in X, system dynamics (cf. Eq. ([7|) 
drive them out to D. Again from a super-macroscopic perspective, we call ^t = dvt{D)/dt the 
outflow of pedestrian from X to D. On the whole, the following conservation principle holds: 

dit dSt ^ /^-.N 

We indicate the rate of emptying of 5, i.e., the rate at which pedestrians appear in X (neglecting 
their spatial distribution), by 

f = /«,/,;«)- 

From the phenomenological point of view, such a rate can depend on several factors such as, for 
instance, the number of pedestrians in X, the amount of agents who still have to enter (i.e., St)^ 
or even the time of the day. 

Coherently, the function / should be designed in order to mimic dynamical arrival processes. 
In the following, an expression for /, suitable to represent smoothly ending queue processes is 
constructed. In particular, we assume: 

• f (X a{St) for some function a : [0, 1] satisfying cr(0) = (i.e., the rate is zero if all 

agents have appeared in X); 



^Notice that intersections of rectangles taken from orthogonal grids, instead, give rise just to rectangular inter- 
sections. 
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• that X cannot be filled up indiscriminately. A certain capacity C > must not be overcome. 
Hence, a logistic factor is considered: / ex (1 — ^). The constant C may satisfy C = pc\L\^ 
pc > being a capacity density threshold of X. In this setting, when If < C 8i mass inflow 
in X from S is allowed, otherwise a reverse flux takes place. 



On the whole, we obtain the expression: 



f{SuIt;t) = a{St 



1-^ 
C 



(22) 



As a matter of fact, a {St) represents the ideal outflow of mass from S to X when It = (i.e., X is 
empty) . If the queue features a constant arrival flux which smoothly ends as N pedestrians enter 
the domain, the function a can be chosen as: 

^St\ St [St 

F being a positive constant and < p < 1 the ratio at which the appearance rate begins to fade 
out. 



cj{St) = F 



L[0,p) 



(23) 



3.4.1 Numerical treatment of the entrance region 



In this section the numerical scheme in Eq. (20) is extended so that pedestrian entering regions 
introduced in Sect. 13.41 can be handled. 



An extended version of domain DUX, including entrance regions, is considered. Therein, a finite 
volume-like partition {Ej}^^-^ consistent witlj^both subsets D and X is introduced. Coherently, 
the discrete density pl^ is defined on such a partition and thus it is defined also on X. 

The dynamic arrival process described by Eqs. (21) and (22) is taken into account in its 
discrete-in-time version, that is: 



Sn-\-l — Sn — Atf{Sn^ ^n+l^ ^n) 

= + ^tf{Sm ^n+i; ^n) — (^n+1 " ^n). 



where J^+i and 4>n+i satisfy respectively 



n+1 



= E 

j:EjCX 

= E 

j-EjCD 



pi+i\Ej\ 



pi+i\Ej\. 



(24) 

(25) 
(26) 



On the whole, the evolution of the crowd through tS, X, and D is computed using the following 
algorithm: 



for n ^ to M do 

Evolvd^p^ on DUX; 



Evaluate <^>n+i and /n+i (using Eqs. (26) and ([25|)); 



Evaluate 5'n+i and i^+i (using Eq. (24)); 
Assigi][^ 



Pn+l 



n+1 - 



end 



Algorithm 1: Procedure to handle pedestrian entering 



^i.e. evolve the discrete densities pn according to algorithm in Eq. ( |20| ) 

^i.e. distribute uniformly the mass /n+i on X - thus overwrite the discrete densities in X. 



li.e. {E, G {E,} 



3Sj = l 



EjCiD^^AEjCiI^^}- 
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Table 1: Parameters used for crowd event simulations 



L B PC N V 

100 m 4 m 1.3 ped/m^ 1500 ped 1.18 m/s 



filling 


?=0.5 


regime 


t= 1.0 


full 


t = 2.0 


walkway 




regime 


? = 3.0 


leaving 


? = 4.0 


regime 


^ = 4.5 



pp 

5 - 



PC 



-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 

x/L 



Pi 






5 1 


15 








0. 


75 1.5 2.25 



filling full walkway leaving 
regime regime regime 



ped /TV 
1 



0.8 



(a) 



(b) 



Figure 8: (a) Instantaneous density fields and recognized regimes - (b) Time history of some crowd 
bulk parameters 



4 Application 

This section is devoted to provide some real world applications of the proposed model. In par- 
ticular, model sensitivity and tunability are discussed in the reference configuration in Fig. [4j 
Moreover, simulation results are discussed for the geometrical the configurations in Fig. [5] 

4.1 Setup overview and some simulated phenomena 

In this section, a prototypical crowd event happening in the straight rectangular walkway D 
depicted in Fig. [4] is considered. The application is defined by assuming that a hypothetical client 
is able to provide some data of the problem setup, particularly L, 5, the expected total number N 



of incoming pedestrians, capacity density pc in Eq. (22), the pedestrian desired speed V — \vd\ 
In the present case values have been selected on the basis of data available in Transportation and 
Civil Engineering literature f5]|TT, 28 and they are summarized in table |4.1[ The length L and 



the desired speed V are reference quantities, thus the time scale T = L/V can be defined as the 
time required to an undisturbed pedestrian to cross the whole facility. In the following, we denote 
by pp the dimensional pedestrian mass density [ped/m?] which is related to the density p through 
the scaling p^ = pN jl? . 

The quantities N and pc are required to characterize the time history of the incoming crowd 



(Sect. 3.4), while the parameters F = 2 and p = 0.9 governing the entrance queue (Eq. 23) have 
been arbitrarily set to simulate a sudden start of the entrance and a smoother end queue, already 
observed in some real world crowd events (e.g. in [4jjTl]). The resulting time history of crowd 
density in the region X is reported in Fig. [sj^b). 

In order to outline the main features of the simulated crowd event. Fig. [sj^a) samples some 
instantaneous density fields and groups them in three recognized regimes: during the filling regime 
pedestrians advance on the partially empty walkway, which is homogeneously filled in the full 
walkway regime. The crowd event gradually ends during the leaving regime. In Fig. [sjb), the time 
history of two integral parameters of the simulated crowd event is also plotted: the number Mf 
of pedestrians along the walkway and the cumulative number Gt of pedestrians walked out, both 
scaled with respect to N. A further bulk parameter can be easily recognized, i.e. the total time 
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of the crowd event defined as 

T„ = inf|t: ^ = l|. (27) 
4.1.1 Sensitivity to free model parameters 

Four free model parameters remains and determine the evolution of the solution, namely the 
constants c, 0, and a. On the one hand, R = 2 m and a = 45°, depending upon the geometry 
of the sensory region, are regarded as case-independent and recoverable from existing literature 
(see e.g., (l0j[28]). On the other hand, the sensitivity of the model to the repulsion constant c and 
to the angle is inquired. The analysis is based on the observable variables T^, which has been 
previously defined in Eq. [27| and 6p: 

Sp= , (28) 

PC 

where ps and pm are the crowd density at the walkway side and at the mid-line, respectively, 
evaluated at the mid-span x = L/2 during the full walkway regime. In other words, such a variable 
is a measure of the chord-wise uniformity of the crowd density, being the span-wise uniformity 
assured during the full walkway regime and the chord-wise symmetry of the solution assured by 
the selected setup. 

In Fig. [9j such variables are plotted versus the dimensionless parameters 2.5 • 10~^ < c < 
12.5 • 10-^, 0° <0 < 5°. 

The crowd event time (Fig. ^s)) is mainly sensitive to the pedestrian-pedestrian repulsion, 
i.e., to c. For the selected incoming pedestrian density pc, Ta approximatively varies from three to 
four times the undisturbed pedestrian crossing time, depending on the value of c. On the contrary, 
Sp (Fig. ^h)) depends upon both parameters and shows both positive values (higher density at 
the walkway sides) and negative ones (high density along the mid-line) in the selected range of 
the model parameters. In other words, for any considered value of T^, there exist a value ^* of the 
angle allowing a nearly homogeneous chord- wise crowd density {6p = 0), while larger or smaller 
values produce nonuniform distributions. In order to detail the chord- wise trend of the crowd 
density and to discuss its phenomenological features, the pp profiles at mid- length are plotted in 



Fig. 10 for different values of and fixed value ofc = 5-10~ . 

For < 6> < 6>*, where ^* ^ 2° with c = 5 • 10~^, the pedestrian- wall repulsion dictated by the 
desired velocity field is lower than the pedestrian-pedestrian repulsion, the chord-wise component 
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Figure 10: Chord- wise crowd density profiles at mid-length for different ^ (c = 5 • 10 ^) 



of the total velocity field is directed toward the lateral walls and crowd density results larger 
at the walkway sides than at mid-chord (Fig. 10, red profiles). Conversely, the pedestrian- wall 
repulsion predominates over the pedestrian-pedestrian one for ^* < ^ < 5°, and the crowd density 
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at mid-chord is larger than at walkway sides (Fig. 10, blue profiles). The crowd density profile 
is almost flat around the value ^* (Fig. [Toj green profiles). In a general modeling perspective: 
i) the parameter d allows one to account for different degrees of repulsion of either both or a single 
walkway side, e.g., a panoramic point on one side only; ii) the parameter c allows one to account 
for different attitudes of pedestrians in accepting the proximity of other walkers, e.g., dictated 
by different travel purposes |5|[28). From a transportation engineering perspective, the parameter 
Q allows one to model the shy distance of pedestrians from the wall |20 and to evaluate the 
effective width of the walkway (e.g., |13|): for 6> = ^* the effective width of the walkway matches 
the geometric width, while shorter effective widths are obtained otherwise (Fig. 
sketches) . 

4.1.2 Tunability of the free parameters 

On the basis of the previous sensitivity analysis a tuning strategy can be suggested. On the one 
hand, the model is compact, in the sense that only two parameters are introduced to characterize 
the two velocity components and to account for distinct repulsion phenomena. On the other hand, 
bulk experimental data directly obtained from real world crowd events are preferable to pointwise 
measurements obtained in laboratory tests, in the perspective of a tuning procedure to be easily 
applied to the engineering practice. 

Hence, we assume the hypothetical client is able to complement the setup data with mea- 
surements obtained on the same walkway of interest or on analogous facilities. In particular, the 
crowd event time Ta is easy to be measured, hence it is usually provided (e.g. in ^^) or estimated. 
Moreover, information about the degree of repulsion of the lateral sides to be adopted is required, 
either qualitatively identified by referring to the prototypical blue-red-green profile shapes in Fig. 
To] or quantitatively expressed by Sp. 

Being such data available, say Ta = 5.2, Sp = 0, the graphs of Fig. [ojcan be used as tuning 
charts: first, the value of the constant c = c is obtained by setting the value of T^/T = in Fig. 
[oja), and then the value of 6 = 6 is recovered by setting the value of dp = dp in Fig. (o^b). In the 
following, the values c ^ 5 • 10~^ and = 0* ^ 2° are retained. 

4.2 Crowd flow along real world footbridges 

The crowd flow along the same real world footbridges described in Sect. |3.2.1[ Fig. [5) is now 
simulated in order to show the ability of the proposed approach to face actual engineering problems. 



The simulations adopt the same setup (Tab. 4.1) and the values of the model parameters set in 



the previous subsection, so that the comparison to the reference test-case (rectangular walkway) 
allows to qualitatively point out the effects of the domain geometry. 
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bottleneck walkway 
(e.g. Chiaves footbridge) 




shifted walkway 

(e.g. Pedro e Ines footbridge) 



curved walkway 

(e.g. Liberty footbridge) 



-0.5 -0.4 -0.3 -0.2 -0.1 



0.1 0.2 0.3 0.4 0.5 

x/L 
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0.75 1.5 2.25 p [ped/m^ ] 



Figure 11: Comparison among real world footbridges: density field during the stationary full 
walking regime 

For instance, Fig. [TT]^a) collects the crowd density instantaneous fields during the stationary 
full walking regime for the considered domains. 

The fields differ significantly both qualitatively and quantitatively: asymmetric patterns arise 
with respect to the chord- wise axis at midspan (bottleneck walkway), to the longitudinal axis 
(curved walkway) or to both (shifted walkway); the maximum crowd density approximatively 
doubles in the bottleneck and shifted walkway respect to the rectangular one. 

Finally, it is worth pointing out that the output of all considered simulations can be equivalently 
interpreted in the statistic way proposed in Sect. |2.3[ In particular, both Figs. [8] and pT] report 
alongside the pedestrian density pp the probability p density associated to the random variables 



5 Conclusions 

In this paper, a model aimed at evaluating the collective evolution of a crowd has been derived 
and analyzed. A model describing the motion of a generic test individual has been firstly obtained 
on a phenomenological basis, incorporating the effect of the perception she has of the collectivity 
in her surroundings. It is worth pointing out that the collective point of view has been treated 
in the general terms of Radon measures, which allow one to distinguish two contributions in the 
perception: an atomic localized one and a continuous distributed one. A procedure has then 
been provided to obtain, out of the equation governing the motion of single individuals, equations 
describing the collective evolution of the crowd. Next, the specific structural class of pedestrian 
footbridges has been considered and the elements required to simulate crowd events therein have 
been given. In particular, a working strategy to construct suitable desired velocity fields, satisfying 
some minimal phenomenological requirements, has been proposed. Issues concerning the equation 
discretization as well as the handling of boundary conditions have been addressed. The reduced 
amount of free parameters (four) present in the final equations, made possible the definition of 
a tuning procedure based on a free parameter sensitivity analysis. To show the ability of the 
proposed approach to face actual engineering problems, crowd events have been simulated in 
different computational domains inspired by real footbridges. 
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